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Abstract. Quantum Hall effect (QHE) devices are a pillar of modern quantum 
electrical metrology. Electrical networks including one or more QHE elements 
can be used as quantum resistance and impedance standards. The analysis of 
these networks allows metrologists to evaluate the effect of the inevitable parasitic 
parameters on their performance as standards. This paper presents a systematic 
analysis of the various circuit models for QHE elements proposed in the literature, 
and the development of a new model. This last model is particularly suited to 
be employed with the analogue electronic circuit simulator SPICE. The SPICE 
macro-model and examples of SPICE simulations, validated by comparison with 
the corresponding analytical solution and/or experimental data, are provided. 
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1. Introduction 

The quantum Hall effect (QHE) has been the basis of resistance metrology for several 
years in DC m and more recently in the AC regime mz]- In DC, networks of several 
QHE elements have been developed to realize quantum Hall array resistance standards 
(QHARS) [51116) ■ intrinsically-referenced voltage dividers [T7| and Wheatstone bridges 
m- In AC, Schurr et al. nni included two QHE elements in a quadrature bridge to 
realize the unit of capacitance. 

In 1988, Ricketts and Kemeny (RK) [5D] proposed a circuit model suitable for the 
symbolic analysis of arbitrary electrical networks containing QHE elements, in DC or 
AC. The predictions of the RK model were also verified experimentally [5Ull21| . Other 
circuit models were developed in later years [211155] and a general method of analysis 
based on the indefinite admittance matrix has been recently proposed m- 

The circuit models proposed in the literature are better suited for the symbolic 
analysis of networks of QHE elements; the matrix method can be employed both for 
symbolic and numerical analyses, typically with the aid of a computer algebra system 
or a numerical computing environment. Nowadays, however, the numerical analysis 
of electrical circuits is commonly carried out by employing analogue electronic circuit 
simulation tools, like SPICE m and its derivatives. SPICE can perform simulations 
of linear and non-linear networks in time and frequency domains, and noise analysis. 
Hence, it could be a very useful tool for the simulation of electrical networks containing 
QHE elements. Unfortunately, even though the above cited circuit models can be 
directly coded in SPICE, several issues arise when trying to run a simulation. The 
aim of this work is to address these issues and to present a working SPICE macro¬ 
model for the simulation of networks of QHE elements. 

In section[2j we review the existing circuit models for QHE elements and present a 
new one. In section [3J we discuss the issues of SPICE modelling and provide a SPICE 
macro-model for an 8-terminal QHE element, which is based on the new circuit model. 
In section H] three simulation examples are worked out in detail: i) a DC analysis 
taking into account parasitic resistances; ii) an AC analysis of a network containing 
a QHE element and a capacitor; hi) a noise analysis. The simulation of example i) is 
compared with an analytical result; those of ii) and iii) are compared with analytical 
and experimental results. [Appendix Aj reports the full derivation of the new circuit 
model. 

2. Brief review of QHE circuit models 

The development of circuit models for the Hall effect is not a recent topic. For models 
related to the classical Hall effect see, e.g., [5TII35| : in this section, instead, we shall 
briefly review the existing circuit models for ideal QHE elements. 

In circuit theory, an n-terminal element is a black-box which can interact with 
its surroundings solely through n conductors called terminals (figure [T|). A basic 
assumption of circuit theory is that an n-terminal element is completely characterized 
by its external behaviour, that is, by the set of admissible terminal voltage-current 
pairs (action at distance is excluded) [361137] . Two elements are thus equivalent from 
the point of view of circuit theory if they have the same external behaviour. 

An electric circuit composed of several interconnected elements can be itself 
turned into an n-terminal element by connecting terminals to n selected nodes of 
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Figure 1. An n-terminal element Af represented with reference polarities and 
directions; terminals are numbered in a clockwise direction; terminal voltages are 
measured with respect to an arbitrary reference point O {datum node)-, terminal 
currents flowing into Af are considered positive. Voltages and currents can be 
functions of time. 


the circuit. If a circuit, viewed as an n-terminal element, is equivalent to a certain 
n-terminal element Af, it is called a circuit model (or equivalent circuit) for Af. 

According to their external behaviour, ideal QHE elements can be classified in 
two categories. In the following, we shall calQ an ideal clockwise (cw) n-terminal QHE 
element with Hall resistance > 0, an element whose terminal voltages and currents 
are related by the equations 

jm — f-m f'm—1 , tn — 1, . . . , n, (1) 

with the convention eg = Cn, and we shall call an ideal eounterclockwise (ccw) n- 
terminal QHE element, an element whose terminal voltages and currents are related 
by the equations 

jm — f-m ^ — 1, . . . , 7f, (2) 

with the convention e„+i = ei. Ideal QHE elements are memoryless, passive (actually, 
they dissipate power) and nonreciproca^ 

The sets of equations (CD and m are, indeed, a consequence of the QHE 
phenomenology, but can also be derived theoretically on the basis of the Landauer- 
Biittiker formalism (see, e.g., [391 Sec. 16.3] and references therein). 

A real QHE device operating at appropriate values of temperature and magnetic 
flux density can be modelled, with a certain approximation, either as an ideal cw or 
an ideal ccw element depending on the orientation of the magnetic field applied to 
the device and on the type of the majority charge carriers. The Hall resistance is 
quantized, i?H = RK/i, where Rk — h/e^ is the von Klitzing constant and i is a 
positive integer called plateau index. In QHE devices realized with GaAs technology, 
the plateau index of interest for resistance metrology is j = 2, so i?H = .Rk/ 2. The 
recommended value of the von Klitzing constant is i?K = 25 812.807443 4(84)12, 

t The terminology has been inspired by the two circuit models shown in the first row of table ITI 
§ They are memoryless because the sets of equations Q and Q relate terminal voltages and currents 
at the same time instant. For definitions of passivity and reciprocity see e.g. IMl Ch. 2]. 
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which yields = Rv^/2 = 12 906.403 7217(42)17 [ID]. The conventional value 
adopted internationally for resistance metrology is i?K -90 = 25 812.8071790, which 
yields i?H -90 = 7?k-9o/2 = 12 906.40351790. In sectionHJ we shall use this last value 
for the Hall resistance of the QHE elements in the example circuits. 

Table [T] reports known circuit models listed in order of appearance in the 
literature: the first row reports the well-known circuit of Ricketts and Kemeny EOj; 
the last one, a new circuit model whose derivation is given in [Appendix A[ For each 
referenced work, the second and third columns show the circuit models corresponding 
to the ideal cw and ccw 4-terminal QHE elements (the cited works might not provide 
both); in order to see how these models can be extended to elements with more than 
4 terminals, the fourth column shows the outlines of the circuit models corresponding 
to a cw 6-terminal element. All circuit models contain only resistors (because of 
memorylessness and passivity) and controlled sources (because of nonreciprocity); 
none of them are related to the actual device physics. 

We remark that the sets of equations (III) and (m are all that is needed to analyse 
circuits containing ideal QHE elements and that the external behaviour of all 
the models listed in table [T] is exactly described by those equations. Nonetheless, a 
circuit model can be a useful basis for taking into account nonidealities (e.g. nonzero 
longitudinal resistances |24l[27|, parasitic elements in AC regime [41]) or to develop a 
SPICE macro-model, as we shall see in the next section. 

3. SPICE modelling 

In SPICE a circuit is represented by a list of statements (netlist), written according to 
a certain specific syntax [DDISD], which describes how the elements that compose the 
circuit are interconnected. A SPICE macro-model (or subcircuit netlist) is a netlist 
with a designated name that can be treated as any other SPICE element to compose 
larger circuits in a hierarchical way. Circuits composed exclusively of resistors and 
controlled sources can be directly translated into SPICE macro-models. Given that 
all the models reported in table [T] are equivalent from the point of view of circuit 
theory, and that many more possessing this property can be conceived, how can we 
select one with the purpose of writing a SPICE macro-model? 

In reality, SPICE modelling requires at least two additional properties, Q and 
m below, but one more, (juil) . would be useful: 

(i) No loops of voltage sources. SPICE requires a circuit that does not contain 
loops composed only of ideal voltage sources (independent or controlled). In fact, 
violating this requirement would lead to two scenarios: either Kirchhoff’s voltage 
law would be violated or the current crossing the loop would be indeterminate 
and the circuit would not have a unique solution. In both cases, SPICE would 
not be able to solve the circuit. 

(ii) No cut sets of current sources. This property is the dual of (|i|). If, in a circuit, 
there is a node where only current sources join, either Kirchhoff’s current law 
would be violated or the node potential would be indeterminate (floating node). 

(iii) Non-dissipative/generative sources. As stated in section [2l QHE elements are 
passive elements that dissipate power. One property that could then be required 
by the circuit model is that all the power is dissipated in the resistors and that 
the overall power absorbed or delivered by the controlled sources is zero. This 
property would also allow SPICE to predict the thermal noise generated by 
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Table 1. List of known circuit models for QHE elements. The first column 
reports references to the original works. The second and third columns report, 
with uniform symbology, the corresponding circuit models for ideal cw and ccw 
4-terminal QHE elements: terminal voltages and currents are chosen according to 
figure^ eij = a — ej; r = Hh/ 2; diamond symbols represent controlled sources. 
In the last column, outlines of the circuit models for a cw 6-terminal element are 
sketched (see the given references for details). 


Reference 


1. Ricketts-Kemeny m 


2. Hartland et al. I22I23I 


3. Jeffery et al. [2111211 


4. Sosso-Capra |25II26| 


5. Schurr et al. 


6. This work 
||Appendix 


4-terminal cw model 4-terminal ccw model 


ei ei 












Outline of the 6-termiiial 
model (cw) 



No generalization has 
been directly proposed by 
the authors 



4 



1 



4 















































Circuit models and SPICE macro-models for quantum Hall effect devices 


6 


the QHE elements correctly (see [Appendix A| for more details on this). Even 
though this property is not strictly required, and even though SPICE has a few 
limitations in noise analysis (e.g. it is not able to analyse correlations directly and 
a workaround is needed |43|i. it might be useful to have a circuit model which is 
as complete as possible. 

By inspection, it is easy to verify that models 1 and 3 of table [T] violate property (P 
and that model 6 violates property ([n]); property P can also be violated by model 5 
if a voltage source is connected between terminals 1 and 2, or 3 and 4. In addition, it 
is straightforward to verify that models 2,4 and 5 violate property (iinlf ||. 

Models 2,4 and 5 cannot be easily modified to satisfy property PU and will not 
be considered further. Models 1 and 3 can be modified to satisfy property P either by 
adding a small series resistance to the loop of voltage sources |33] or by removing one of 
the controlled sources in the loop (this does not alter the model’s behaviour). However, 
none of the above two remedies is completely satisfactory: the former because the effect 
of the additional resistance on the simulation accuracy cannot be easily predicted in 
the case of many interconnected QHE elements; the latter because it destroys the 
symmetry of the model, making it more difficult to develop possible refinements to 
take account of nonidealities. Eor model 6, it is worth noting that the potential e 
of the floating node can be fixed arbitrarily without altering the model’s behaviour; 
this allows the satisfaction of property p]) without the need for additional elements 
and without destroying the circuit symmetry. Therefore, model 6 was selected to 
implement a SPICE macro-model. 

Listings d] and H] report, respectively, the SPICE macro-models for ideal cw and 
ccw 8-terminal QHE element^. These macro-models are a direct translation of the 
general circuit models derived in [Appendix A[ The actual number of terminals is nine: 
the ninth terminal, labelled C, corresponds to the node joining the controlled current 
sources and has to be connected to an arbitrary potential (e.g. ground). The default 
value of i?H is 1 H, but other values can be assigned when calling the macro-model 
(see next section for an example). 

Listing 1. SPICE macro-model for an ideal cw 8-terminal element. 

.subckt qheScw 12345678C params: RH=1 
R1 1 2 ■{2*RH} 

R2 2 3 ■{2*RH} 

R3 3 4 {2*RH} 

R4 4 5 {2*RH} 

R5 5 6 {2*RH} 

R6 6 7 {2*RH} 

R7 7 8 {2*RH} 

R8 8 1 {2*RH} 

G1 1 C 2 8 {1/(2*RH)} 

G2 2 C 3 1 {1/(2*RH)} 

G3 3 C 4 2 {1/(2*RH)} 

G4 4 C 5 3 {1/(2*RH)} 

G5 5 C 6 4 {1/(2*RH)} 

II For instance, for model 5 (cw), the power delivered by the controlled sources is p = ‘2r{jij2 -hjiji), 
which is generally not zero. 

1 These macro-models should work out of the box, or with just slight modifications, with any modern 
SPICE simulator that accepts subcircuits with parameters. In this work, we have chiefly used LTspice 
IV from Linear Technology Corporation but tests have been carried out also with the open source 
simulator Ngspice |46L with PSpice A/D from Cadence Design Systems and with TINA-TI from Texas 
Instruments. The authors can provide advice on adapting the circuit models and the simulations here 
described to other SPICE-based analogue circuit simulators. 
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G6 6 C 7 5 {1/(2*RH)} 
G7 7 C 8 6 {1/(2*RH)} 
G8 8 C 1 7 {1/(2*RH)} 
. ends 


Listing 2. SPICE macro-model for an ideal ccw 8-terminal element. 


.subckt 

qhe8ccw 12345678C params 

R1 1 

2 

■{2*RH} 

R2 2 

3 

■{2*RH} 

R3 3 

4 

{2*RH} 

R4 4 

5 

{2*RH} 

R5 5 

6 

{2*RH} 

R6 6 

7 

{2*RH} 

R7 7 

8 

{2*RH} 

R8 8 

1 

■{2*RH} 

G1 C 

1 

2 8 {1/(2*RH)} 

G2 C 

2 

3 1 {1/(2*RH)} 

G3 C 

3 

4 2 {1/(2*RH)} 

G4 C 

4 

5 3 {1/(2*RH)} 

G5 C 

5 

6 4 {1/(2*RH)} 

G6 C 

6 

7 5 {1/(2*RH)} 

G7 C 

7 

8 6 {1/(2*RH)} 

G8 C 

8 

1 7 {1/(2*RH)} 

. ends 



4. Examples 


f.l. Double-series interconnection of two QHE elements 


Several QHE elements can be interconnected to obtain multiples, submultiples 
or fractions of the Hall resistance. Multiple-series, parallel m and bridge 
connections |48] can be employed to reject the effect of the inevitable contact and 
wiring resistances. In this section, SPICE is used to analyse, in the DC regime, the 
effect of the parasitic resistances in a double-series connection; the result of the SPICE 
analysis is then compared with an analytical solution obtained with the technique 
described in m- 

Figured] shows the complete circuit diagram for the analysis of the double-series 
connection. For ease of analysis, only two parasitic resistances, ri = eir and r 2 = e 2 r 
(r = Ryi/2), are taken into account. The resistance i?g of the double-series is 
calculated between the terminal la and ground, i?g ' = Vi^/I q-, where Iq can be 
chosen of 1A to simplify the calculation. The effect of the parasitic resistances can be 
quantified by the relative discrepancy 


2i?H 


(3) 


where 2i?H is the resistance of the double-series when ri = r 2 = OH. At the third 
order in ei and £ 2 , the analytical solution yield^ 

rtheo _ ^1^2 _ efe2 + eie 2 , . 

16 64 ■ ^ ’ 


■*" The authors can provide the Mathematica® notebook used to obtain the analytical solution. 
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Figure 2. Circuit diagram for the DC analysis of a double-series connection. 
r\ and r 2 represent the wiring and contact resistances between the two devices; 
other contact resistances are neglected for simplicity. 


Listing 3. SPICE netlist corresponding to the circuit of figure [2] The included 
file qheScw.sub should contain the SPICE macro-model of listing^ In this 
example we have set ei = 0.15t and £2 = 0.35t, where t is a dimensionless scaling 
parameter. The dot command .op declares that SPICE should perform a DC 
operating point analysis. The dot command .meas (LTspice IV specific) performs 
the calculation of the parameter (5 as defined by 0. 

QHE double - series circuit simulation 

* Includes the macro-model 
.inc qheScw.sub 

* Definition of circuit parameters 
.param RH=12906.4035 r={RH/2} 

.param t =0.01 

.param r1={0.15*t*r} r2={0.35*t*r} 

* Circuit netlist 

XUl lA 2A lA 4A 5A 6A 7A 8A 0 qheScw 
+params: RH={RH} 

XU2 IB 2B 3B 4B 0 6B 0 SB 0 qheScw 
+params: RH={RH} 

10 0 lA 1 
rl 5A IB {rl} 
r2 7A 3B {r2} 

* Analysis directives 
. op 

.meas op delta param (V(1A)-2*RH)/(2*RH) 

. end 

The SPICE netlislQ which corresponds to the circuit of figure [5] is reported in 

* Modern SPICE simulators come with handy graphical schematic editors which alleviate the chore 
of writing netlists. Here, however, we have preferred to provide netlists, using basic elements and 
directives, which can be easily converted from one SPICE dialect to another (the dialect employed 
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Table 2. Comparison, for the circuit of figure[2]with ri = 0.15tr and r 2 = 0.35t r, 
between the relative discrepancy obtained from ((Sjl, and <5®'™, obtained by 

simulating with SPICE the corresponding netlist of listing |3] The third column 
^sim ].gpo]-^,g results obtained with the LTspice IV’s default solver; similar 
results are obtained with other SPICE-based simulators. The column 5®™ (alt) 
reports the results obtained with the LTspice IV’s alternate solver with reduced 
round-off errors. 


t 

jtheo 


(alt) 

10-1 

3.24 X 10-3 

3.24 X 10-3 

3.24 X 10-3 

10-2 

3.28 X 10-1' 

3.28 X 10-1' 

3.28 X 10-1' 

10-3 

3.28 X 10-® 

3.28 X 10-® 

3.28 X 10-° 

10-^ 

3.28 X 10-11 

5.62 X 10-11 

3.28 X 10-11 

10-3 

3.28 X 10-13 

-1.70 X 10-1° 

3.22 X 10-13 


listing [31 with example parasitic resistances ri = 0.15tr and r 2 = 0.35tr, where we 
have introduced a scaling parameter t to check the round-off errors due to floating 
point arithmetic and numerical algorithms. From ([31), i5 is expected to scale as f^. 

Table H] reports a comparison, for different values of t, between the values of 
^theo obtained from (|31) and those of obtained from the SPICE simulation of 
listing 131 The third column of table [31 can be considered representative of the 
results obtainable with common SPICE simulators: as can be seen from the reported 
results, when S falls below 10“^°, diverges from The last column 

(alt) reports, instead, the results obtained with the LTspice IV’s alternate solver, 
which, according to the LTspice IV manual, has reduced round-off errors. In any 
case, whatever the simulator employed in this kind of analysis, we suggest to write all 
parasitic resistances as functions of a scaling parameter to keep under control round-off 
errors. 

4-2. QHE gyrator 

QHE elements, like their classical counterparts, can be employed as gyrators |15ll50| . 
At the angular frequency w, the circuit of figure[31 which includes the load capacitor Cl 
with negative reactance Vp = — (wCl)~^, realizes between terminal I and ground an 
impedance Z(uj) = Vijlo = i?(cLi)-|-jV(a;) having positive reactance, that is, X{uj) > 0 
for any w. Therefore, at any given frequency, the circuit of figure [3| behaves like an 
RL two-terminal element (the values of R and L are frequency-dependent). 

Network analysis yields 

Z{io) = i + ^clRl + iu;C^Rn] , (5) 

for which X{uj) = Im Z{u!) > 0. 

The behaviour of such an unconventional QHE circuit was checked experimentally 
with an 8-terminal GaAs-AlGaAs Hall bar working at the temperature T = 1.6 K and 
plateau index i = 2. Terminals 3 and 7 were connected to a variable capacitance 
box; terminals pairs (1,2) and (5,6) were connected in double-series configuration to 
an LCR meter (Agilent mod. 4284A) which performed the impedance measurements. 


here is mainly that of LTspice IV). 
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Figure 3. Circuit diagram for the AC analysis of a QHE element used as gyrator. 
The impedance of interest is Z{u)) = Vi/Iq. 



Figure 4. Circuit diagram for the AC analysis of a QHE element used as gyrator 
in the experimental set-up described in the text. Cpi and C'p 2 represent the cable 
capacitances. The LCR bridge measures the impedance = Vi/Iq^ which is 

different from the impedance Z{uj) defined in figure [3| 


The measurement frequency / = a;/(27r) was chosen at 1233 Hz, so that the reactance 
of a 10 nF capacitor is Xl ~ —Ra- 

However, the above described experimental set-up is not accurately modelled by 
the circuit of figure [3] and by the corresponding equation ([5]), because in figure |3] 
the cable capacitances and the real operation of the LCR meter are not taken into 
account (see [7j for an example about the effect of the cable capacitances in AC 
QHE measurements). The experimental set-up is instead better modelled by the 
circuit of figured where Cpi and Cp 2 represent the cable capacitances, and where the 
transimpedance = Vi/Iq represents the impedance actually measured by the 

LCR meter, that is, the ratio of the meter’s high-side voltage to the meter’s low-side 
current. 

Even though the circuit of figure 3] can be solved analytically, this kind of analysis 
is not straightforward, and SPICE simulation can provide a quicker response. Listing!?] 
reports the SPICE netlist for the AC analysis of the circuit of figure (H for different 
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values of Cl at the fixed frequency of 1233 Hz. The given values for Cl are those used 
in the experiment: OnF, InF, 3nF, 7 nF, 10 nF, 15 nF, 20 nF and 30 nF. The values 
of the parasitic capacitances reported in the listing, Cpi = Cp 2 = 300 pF, are just 
rough estimates and not the result of a measurement. 

Listing 4. SPICE netlist corresponding to the circuit of figure [3| The . step 
directive is used to automatically run simulations with different values of Cl ; the 
.ac directive declares that SPICE should perform an AC sweep analysis, here of 
just one frequency point; the .probe directive is used to save node voltages and 
branch currents for further analysis, plots etc. 

QHE gyrator simulation 

* Includes the macro-model 
.inc qheScw.sub 

* Definition of circuit parameters 
.param RH=12906.4035 

.step param CL list 0 In 3n 7n lOn 15n 20n 30n 

* Circuit netlist 

XUl 113455780 qheBcw 
+ params : RH = {RH} 

10 0 1 AC 1 
CL 3 7 {CL} 

Cpl 7 0 300p 
Cp2 3 0 300p 
VG 5 0 0 

* Analysis directives 
.ac lin 1 1233 1233 
.probe V(l) I(VG) 

. end 

Results are reported in figure^ where the measured values of -R(C'l) and V(C'l) 
at the frequency of 1233 Hz are compared with the values obtained from the SPICE 
simulation of listing |4] and with those obtained from ([Q • The agreement between the 
measured values and those obtained from the SPICE simulation is within 1 %. Further 
model refinements can be easily implemented in the SPICE simulation. 

4.3. Noise analysis 

Noise is a fundamental limit to the resolution of a measuring system. Noise analysis 
is therefore a useful tool to predict the resolution achievable in a measurement. 
SPICE has built-in noise models for resistors and semiconductor devices which allow 
a complete noise analysis of circuits, taking into account all the main noise sources. 
In the case of circuits containing QHE elements, the noise analysis can include 
the thermal noise generated by these elements and the noise generated by possible 
auxiliary elements, such as voltage or current sources, detectors etc. 

In this section we present the noise analysis of the circuit of figure [SI for which the 
thermal noise properties have been investigated both theoretically and experimentally 
in m- The one-sided cross-spectral density function Sah{f) between the voltages Va 
and Vh due to thermal noise is S'ab(/) = 2A:BTi?H [SI]j where fee is the Boltzmann’s 
constant and T is the thermodynamic temperature of the QHE element. 

Figure |7| shows the circuit needed to evaluate <S'ab(/) with SPICE. Since SPICE 
cannot evaluate directly cross-spectral density functions, but only (auto-) spectral 
density functions, a trick should be employed [53]. Four voltage controlled current 
sources, with transconductances of 1 S, generate two currents, one proportional to the 


R, X/kn 
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Figure 5. QHE gyrator: comparison of results from experiment (+), simulation 
(o) and equation Q (solid line), for different values of the load capacitance Cl 
at the fixed frequency of 1233 Hz. 



Figure 6. Example circuit for noise analysis. The objective is to evaluate the 
cross-spectral density function S'ab(/) between the voltages Ua and Ub- 


sum of Ua and the other proportional to their difference; it can be shown that 

5ab(/) = t[5+(/)-5_(/)], (6) 

where S'+(/) is the spectral density function of Ua+Ub and S-{f) is the spectral density 
function of Ua — Wb- Both S+{f) and S-{f) can be directly evaluated by SPICE. The 
additional, noiseless, voltage source Vq is actually ineffective in this analysis, but it is 
required by SPICE to have a fictitious input. 
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Figure 7. The circuit of figure [6| with the additional elements needed for the 
SPICE evaluation of S'ab(/)* 


Listing [S] reports the SPICE netlist for the noise analysis described above. The 
subcircuit probe_correlation is used to calculate the sum and the difference of Ua 
and in two successive steps: when the parameter gb is set to 1 the output node 
(out) of probe_correlation yields the sum of the two voltages; when gb is set to 
-1, the difference is instead obtained. The simulation temperature is set to 2K. The 
calculation of 5’ab(/) or of the ratio 5'ab(/)/(4A:Bi?H^) can be done automatically by 
the SPICE graphical post-processor, but we shall not dwell on the details here. The 
result obtained from the simulation coincides with the theoretical one. 

Listing 5. SPICE netlist corresponding to the circuit of figure [3 
QHE noise simulation 

* Includes the macro-model 
.inc qheScw.sub 

* Defines a subcircuit for calculating 

* the sum and the difference of two 

* voltages 

.subckt probe_correlation A B OUT 
G3 0 J A 0 1 

G4 0 J B 0 {gb} 

V2 0 J 0 

H2 OUT 0 V2 1 

.ends probe_correlation 

* Definition of circuit parameters 
.param RH=12906.4035 

.temp -271.15 

.step param gb list -1 1 

* Circuit netlist 

XUl 123406700 qheScw 
+ params : RH = {RH} 
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VO 1 0 0 

XXI 2 4 OUT probe_correlatioii 

* Analysis directives 

.noise V(OUT) VO dec 10 100 10k 


. end 


5. Conclusions 

In this paper we reviewed the modelling of the ideal quantum Hall effect element as a 
n-terminal electrical network. The unique set of network equations can be represented 
with different equivalent circuit models published in literature. An analysis of these 
models shows that they are unsuitable to be directly coded in circuit simulation 
software such as SPICE, because they generate singular representations or because 
of errors in the noise simulation. A new QHE circuit model, specifically designed 
to be implemented in SPICE, is proposed. This paper provides several examples of 
electrical circuits including QHE elements, of their coding in SPICE and hints for 
proper simulation runs. For each case, the simulation outcome is compared with the 
results from analytical modelling and also, in two cases, with experimental data. 
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Appendix A. Yet another model 
Appendix A.l. Model derivation 

In this section we derive model 6 of table [T] for an ideal n-terminal QHE element on 
the basis of the analysis method described in [5S] and briefly reviewed below. In the 
derivation, sinusoidal regime is assumed: voltages and currents are to be understood 
as voltage and current phasors and are denoted with capital letters. Labelling of 
terminals and reference directions are again that of figure [TJ 

The sets of equations © and m can be rewritten in matrix form as |29| 


(A.l) 


J = Y,E, 


where J = {Ji,..., Jn)"^ (the superscript T denotes transposition) and E = 
{El, ..., En)^ are column vectors, and where the n x n indefinite admittance matrix 
(see for a definition) Y, is either 


/I 0 0 -I \ 

-I I 0 ••• 0 


1,CW — 


I 

i?H 


(A.2) 


0 -I I 0 
V 0 0 -I I / 
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in the case of an ideal cw n-terminal element, or 


/ 1 -1 


Y — 

^ 1,CCW - 


1 

i?H 


0 


1 -1 


0 \ 
0 


(A.3) 


0 ••• 0 1 -1 

V -1 0 ••• 0 1 y 

in the case of an ideal ccw n-terminal element. 

The two above matrices are real matrices. A basic theorem of algebra |541 Ch. 1] 
states that every matrix can be decomposed into the sum of a symmetric part and an 
antisymmetric one: 


Yi = Yf 


Yf 


with 


Y^ = -{Yi + YY) and 




(A.4) 


(A.5) 


Since, by definition, {Y^)ki = {Y^)ik, the matrix is associated to a reciprocal 
network [551 Ch. 16]; in addition, this associated network can be realized with just 
resistors because 5]® is also real. The antisymmetric part Y^^ is instead associated to 
a nonreciprocal network that can be realized with voltage controlled current sources. 
For the matrices in (|A.2I) and (|A.3ll the symmetric and antisymmetric parts are 


y:® = y:® 


1 

i?H 


/ 1 - 


V -5 
/ 0 


Yt, = 


1 

i?H 


V I 




(A.6) 


1 y 


(A.7) 


0 y 


and 


y-A _ _y;A 
i,ccw i,cw 


/ 0 - 


1 

i?H 


0 




0 


(A.8) 


0 y 


From taking into account (IA.4|) . we obtain that the terminal currents can 

be decomposed into the sum of two contributions, one associated to 1^® and one to 
KA; 


J = J® 


(A.9) 














Circuit models and SPICE macro-models for quantum Hall effect devices 


16 


Em 



Figure Al. Portion of the circuit model associated with the symmetric part 
of Yi. 


with 


jS = and . (A.IO) 

This means that and can be considered associated to two electrical networks 
connected in parallel. 

Let J® = (jf,..., J®)"^: from (lA.lOl) and (|A.6|) the cnrrent at the mth terminal 
is (m = I,... ,n-, Eo = En and A„+i = Ei) 


jS _ Em — Em-lf‘2 + Am+l/S 
-Kh 

_ Em Em—1 ^ Em Em+1 


2i?H 


2i?H 


(A.ll) 

(A.12) 


The last equation above corresponds to a network where each terminal is connected to 
its adjacent terminals through a resistor with resistance 2i?H, as shown in figure lAT] 
This can also be obtained directly from the properties of the indefinite admittance 
matrices [5Sl Ch. 2, Sec. 2.2]. 

Now, let = {Jij ..., from (IA.10|) . (IA.7I) and l|A.8p . the current at the 

mth terminal is 


_ I Em+1 Em—1 


(A.13) 


where the plus sign holds for a cw element and the minus for a ccw one. This 
equation can be obtained with a voltage controlled current source, driven by the 
voltage difference Em+i — Em-i, which draws from the mth terminal the current 
±{Em+i — Em-i)/{‘IRu) injecting it into an additional internal node E (figure|A2|. 

Finally, combining in parallel the two networks, one obtains the circuit model 
of figure IA3I This circuit model obviously satisfies property Q ; property (|n]) can be 
satisfied by fixing the potential E of the central node of figure IA3I to an arbitrary 
value. 


Appendix A.2. Power 


To prove that the circuit model of figure [AH] satisfies property (Iml) . let us recall that, in 
sinusoidal regime, the average power entering an n-terminal element is given by [551155] 





(A.14) 
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Em 



Figure A2. Portion of the circuit model associated with the antisymmetric part 
of Vi; Eij = Ei — Ej-y for the voltage controlled current sources, the plus 
sign holds for a cw element and the minus for a ccw one; E is the potential of an 
additional internal node. 


Em 



Figure A3. Portion of the final circuit model obtained by combining in parallel 
the two circuits of figures ED and EH 


= l(J*E + E*J), (A.15) 

where the operator Re denotes the real part and the asterisk * denotes the conjugate 
transpose. Substituting (lA.ll) in (|A.15I) . and taking into account that is real, yields 


^E*{Y*+\-)E, 

(A.16) 

^E*{Y;^ + ^)E, 

(A.17) 

-E*Y^E. 

2 

(A.18) 


Equation (IA.18I1 says that the power dissipated in the ideal QHE element is related 
only to the symmetric part of the corresponding indefinite admittance matrix. Thus, 
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no power is delivered or absorbed by the circuit associated with the antisymmetric 
part, and property ((mil is satisfied. 

Appendix A.3. Noise modelling 

When SPICE performs the noise analysis of a circuit, it automatically assigns noise 
sources to certain circuit elements, in accordance with known noise models. The 
insertion of arbitrary noise sources is not permitted. In the case of the macro-models 
of listings [T] and [21 SPICE will assign noise sources to the resistors only: the controlled 
sources are considered noiseless. To each resistor, SPICE assign a noise source which 
generates the corresponding thermal (equilibrium) noise at the simulation temperature 
(specified by the SPICE parameter TNOM, 27°C by default). Therefore, to prove 
that the satisfaction of property (Iml) allows SPICE to predict the equilibrium noise 
generated by the QHE elements correctly, we need only to prove that the equilibrium 
noise of a QHE element coincides with that of the resistive network associated with 

The equilibrium noise properties of a general passive multiterminal network 
were derived by Twiss m from thermodynamic considerations. Later, Biittiker 
derived an equivalent result for multiterminal conductors from a quantum scattering 
theory [liiniji. These results were verified experimentally by the authors in [51) . 

For the present purposes, the result given by Twiss is in a more convenient 
form. This result can be stated as follows m- consider a passive n-terminal 
element and assume that all the terminals are grounded and one, say terminal n 
without loss of generality, is considered as reference; then, the one-sided cross-spectral 
density function m Sj^j^{f) between the noise currents jp and jq at the terminals 
p = 1,..., n — 1 and ( 3 '=l,...,n — lis given by 

if) = ‘^kBTiYpq + Y*^) , (A.19) 

where fee is the Boltzmann’s constant, T is the thermodynamic temperature and 
Y = (Ypq) is the short-circuit {n — 1) x (n — 1) admittance matrix obtained by 
considering the n-terminal element as an (n — l)-port with ports defined between 
the terminals 1,..., n — 1 and the reference node n. 

The short-circuit admittance matrix can be obtained by the indefinite admittance 
matrix associated to the n-terminal element by deleting the nth row and column ESI; 
then, for and ideal QHE element, Ypq = (Y{)pq and 


2kBT{Ypq + Yq*p), 

(A.20) 

2k^T [{Y.)pq + (y)*p] , 

(A.21) 

2kBT [{Y.)pq + i\\)qp] , 

(A.22) 

4kBTiY^)pq, 

(A.23) 

2kBTUY%q + iY^)qp . 

(A.24) 


Hence, and taking into account that the reference node is arbitrary, the equilibrium 
noise of a QHE element coincides with that of the resistive network associated with 


The two approaches are complementary: the derivation of Twiss can be considered analogous to 
that of Nyquist m for the Johnson-Nyquist noise generated by a resistor; the derivation of Biittiker 
can be considered analogous to that of Landauer 1611 . 
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